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Abstract 

We study the structure of antiproton spectra at extreme subthreshold bom- 
barding energies using a thermodynamic picture. Antiproton production pro- 
cesses and final state interactions are discussed in detail in order to find out 
what can be learned about these processes from the observed spectra. 
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I. INTRODUCTION 



There exist numerous examples for the production of particles in heavy ion collisions at 
bombarding energies well below the single nucleon-nucleon threshold [Q. This phenomenon 
indicates collective interactions among the many participating nucleons and thus is expected 
to give information about the hot and dense matter formed in these collisions. At beam 
energies around 1 GeV per nucleon the most extreme of these subthreshold particles is the 
antiproton. Therefore, it is believed to be a very sensitive probe to collective behaviour in 
nucleus-nucleus collisions. 

However, presently neither the production mechanism nor the final state interactions of 
antiprotons in dense nucleonic matter are well understood. The antiproton yields measured 
at GSI and BEVALAC seem to be described equally well by various microscopic models 
using different assumptions about the production mechanism and particle properties in dense 
nuclear matter This ambiguity raises the question which kind of information can be 

really deduced from subthreshold p spectra. In this paper we use a simple thermodynamic 
framework as a background on wich we can systematically study the influence of different 
assumptions on the final p spectrum. 

In the next section we will focus on the production process. Following a discussion of 
the final state interactions of the antiproton in dense hadronic matter in Section |T|, we use 



in Section a one-dimensional hydrodynamic model for the exploding fireball to clarify 
which features of the production and reabsorption mechanisms should survive in the final 
spectra in a dynamical environment. We summarize our results in Section 0. 

II. PRODUCTION OF ANTIPROTONS IN HEAVY ION COLLISIONS 

A. The Antiproton Production Rate 

Unfortunately very little is known about the production mechanism for antiprotons in 
dense nuclear matter. Therefore, we are forced to use intuitive arguments to obtain a 
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plausible expression for the production rate. As commonly done in microscopic models |§] 
we consider only two body collisions and take the experimentally measured cross sections for 
p production in free NN collisions as input. The problem can then be split into two parts: 
the distribution of the two colliding nucleons in momentum space and the elementary cross 
sections for antiproton production. The procedure is later generalized to collisions among 
other types of particles (Section [II C|) using phase space arguments. Formally, the antiproton 
production rate P, i.e. the number of antiprotons produced in the space-time cell d'^x and 
momentum space element d^p, is given by 

P = T.J ds2w{s)^^{s) j du,du,5{s-{p.+p,f)Uf,. (1) 

In this expression the summation includes all incoming particle species i and j with 4- 
momenta Pi and Pj and phase space distributions fi and respectively. The integration is 
over the invariant collision energy s = {pi + Pj)"^ and the phase space of the two colliding 
particles: 

du;,^-^^6{p^~m^)Q{p^). (2) 

Finally, w{s,mi,mj) = \J {s — {rrii — rrijY) (s — (mj + rrij)'^) is the well known flux factor, 
and d^aij^p{s)/d^p the differential p production cross section. The influence of the shape of 
the phase space distributions and of the differential production cross sections will now be 
investigated in more detail, beginning with the case of p production by A^A^ collisions. 

B. Phase Space Distributions and Parametrisation of the Differential Production 

Cross Sections 

Assuming that the momentum distribution of the nucleons gradually evolves from 5-peaks 
at the momenta of the colliding nuclei to a relativistic Maxwell-Boltzmann distribution, an 
analytical expression for the phase space distribution fij was derived by Schiirmann et al. in 
Ref. |T^. We use this model as the basis for a numerical evaluation of the probability A(s, t) 
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of finding two nucleons at a center-of-mass energy ^/s at time t (the second integral on the 
right hand side of Eq. (|l])). Fig. |l| shows A(s,t) at t = 0.2, 2, 5 fm/c after the nuclei have 
started overlapping, assuming a bombarding energy of 1 GeV/A. Although the initial 5-like 
distribution at s = 5.5 GeV^ widens quite rapidly, the extreme states needed for antiproton 
production [s > (4m)^ = 14.1GeV^) are populated substantially only after about 5 fm/c. At 
this moment the distribution is already nearly a relativistic Maxwell-Boltzmann distribution 
(dashed line) which maximizes the p production rate in our approach. This fact motivates us 
to neglect pre-equilibrium production and to switch on p-production only when the system 
has reached local thermal equilibrium. Similiar conclusions have been drawn from the results 
of a microscopic kinetic model for the phase space evolution |Q . There it was demonstrated 
that only nucleons which have suffered at least several collisions are effective in p production. 
There are no experimental data for the nucleonic processes 

p+p^p+p+p+p 



n + p^n + p + p + p 



close to the production threshold. As already pointed out in [^, different extrapolations 
of the data at higher energies down to threshold change the p-yield by orders of magnitude. 
However, we will show that not only the total yield, but also the shape of the spectrum is 
very sensitive to the threshold behaviour. We choose for the production cross section the 
parametrization 

(v/i-4m)" 

cr{s) = ctq — , (3) 

where the parameter a models the shape of the cross section near threshold. The normal- 
ization (To is determined by the experimentally measured point closest to threshold. If the 
behaviour of a{s) were dominated by the available phase space for the outgoing particles 
(this is the fundamental assumption in Fermi's statistical model [T^), then a = 7/2. 
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In order to obtain from the total cross section (j^) a formula for the differential cross 
section, we assume like others that the momentum distribution of the produced particles is 
mainly governed by phase space. This leads to the simple relationship 

Here Rn is the volume the n-particle phase space, which can be given analytically in the 



non-relativistic limit ||T^. P is the total 4- momentum of the 4-particle final state, which 
reduces to (i/i, 0) in the CMS of the four particles. 

Inserting (|^) and @ into (|l]) and using thermal phase space distributions, the production 
spectrum can be evaluated numerically. As seen in Fig. the sensitivity on the threshold 
power a is particularly pronounced for low-energy antiprotons. The strong deviation at low 
p energies from a thermal shape (dashed lines) for small values of a is most remarkable. 

Using the Fermi value a = 7/2 for the cross section, an analytical expression for the 
production rate can be obtained: 



Here the densities of the incoming particles, YSmin(p) is the threshold for produc- 

ing an antiproton with momentum p, and T is the local temperature of the environment. 
Clearly the rate P is proportional to the product of the densities of the two colliding par- 
ticle species i and j. More importantly, this formula shows the extreme sensitivity on the 
production threshold and on the medium temperature which both appear in the exponent 
of the exponential function. 

C. Production by Resonances and Mesons 

The production of resonances as intermediate energy storage has been pointed out as an 



important mechanism for p production in |]T^ . Because no experimental data are available for 



antiproton production by baryonic resonances, we follow common strategy and assume 
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that the cross section aij^NNNp{s) is independent of the internal state of excitation of the 
colhding baryons in our thermal picture. The consequences of this assumption are quite 
obvious. While the distance to the p-threshold is reduced by the larger rest mass of the 
resonances, the mean velocity of a heavy resonance state in a thermal system is smaller 
than that of a nucleon. Both factors counteract each other, and indeed we found that the 
total rate P is not strongly changed by the inclusion of resonances. 

The role of pionic intermediate states for p-production in pp-coUisions was pointed out 



by Feldman [|T^]. As mesons are created numerously in the course of a heavy ion collision, 
mesonic states gain even more importance in this case. In fact, Ko and Ge ||13| claimed that 
PP PP should be the dominant production channel. Relating the pp-production channel 
to the PP annihilation channel [jl3[ by 

where 5* = 1 is the spin factor for the p, the production rate can be calculated straightfor- 
wardly from Eq. (|ip: 

gp TtT 16 2 UT72^ r. / 



The spin-isospin factor of the p is gp = 9, and E is the energy of the produced antiproton. 
The modified Bessel function Ki results from the assumption of local thermal equilibration 
for the p-distribution. Expanding the Bessel function for large values of 2E/T we see that 
the "temperature" Tp of the p-spectrum is only half the medium temperature: Tp = |T. 
Finally, the meson-baryon collision process 

mi + Bj ^ Bk+p + p (7) 

remains as the last prototype of collisions with two incoming particles. We choose a 



parametrisation for the cross section which reflects its qualitative behaviour WM: 



Comparing this form with measured data on vr j9 — npp colhsions [O], a value of a. 



0.35 mb is obtained. Due to the threshold behaviour of Eq. (W) and the rather large value 



of a^j, it turns out [|I^ that this process is by far the most important one in a chemically 



equilibrated system. However, this chemical equilibration - if achieved at all - is reached 
only in the final stages of the heavy ion collision when cooling has already started. So it is 
by no means clear whether the dominance of the meson-baryon channel remains valid in a 



realistic collision scenario. This point will be further discussed in Section |V 



III. FINAL STATE INTERACTION OF THE ANTIPROTON 

Once an antiproton is created in the hot and dense hadronic medium, its state will be 
modified by interactions with the surrounding particles. Two fundamentally different cases 
have to be distinguished: elastic scattering, which leads to a reconfiguration in phase space, 
driving the momentum distribution towards a thermal one with the temperature of the 
surrounding medium, and annihilation. Each process will be considered in turn. 

A. Elastic Scattering 

The time evolution of the distribution function f{p,t) is generally described by the 



equation |18 



/(p,t2) = j w{p,p'-MM)f{p'M)dp' (9) 

where p'; ^2, ii) is the transition probability from momentum state p' at time ti to state p 
at t2. Because the number density of antiprotons is negligible compared to the total particle 
density in the system, the evolution of /(p, t) can be viewed as a Markov process. Assuming 
furthermore that the duration of a single scattering process r and the mean free path A are 
small compared to the typical time scale 5t and length scale 5r which measure the variation 
of the thermodynamic properties of the system. 



Eq. can be transformed into a master equation. Considering the structure of the differ- 
ential pp cross section one notices that in the interesting energy range it is strongly peaked 



in the forward direction |p!9|-[2T|]. Therefore, the master equation can be approximated by a 
Fokker-Planck equation p2| : 

(10) 



For the evaluation of the friction coefficient A and the diffusion coefficient D we follow the 
treatment described by Svetitsky |]23|. For the differential cross section we took a form 
suggested in W^: 



da 
It 



{h + cA) 





(a + A)^ 


t\ 




t 





(11) 



where the constants a, 6, c are fitted to experimental data. A is the wave number of the 
relative motion, t the transferred four momentum and Ji a spherical Bessel function. The es- 
sential result [jl^ is that the coefficients A and B are nearly independent of the p momentum 
and satisfactorily fulfill the Einstein relation 



D = mTA. 



(12) 



The effect of these coefficients can be visualized by examining the time evolution of a 
non-relativistic Maxwell-Boltzmann distribution whose initial slope Tq is different from the 
temperature of the medium: 



= 0) =iVe-P'/2mTo_ 



(13) 



In a medium at fixed temperature T = 100 MeV and normal nuclear density po the solution 
f{p,t) of the Fokker-Planck equation can be calculated analytically: 

^2 At 



f{p,t) ~ exp 



■p 



exp 



V/2mTeff(t) 



(14) 



2D /A (e2^* - l) + 2mTo 

This shows that the exponential shape of the distribution function is maintained throughout 
the time evolution, but that the slope Tefj(t) gradually evolves from Tq to the value D/mA 
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which, according to the Einstein relation ([T2D, is the medium temperature T. Looking at 



Fig. 1^ it is clear that after about 10 fm/c the spectrum is practically thermalized. Therefore 
initial structures of the production spectrum (like the ones seen in Fig. |^) are washed out 
quite rapidely, and their experimental observation will be very difficult. 

B. Annihilation 

The annihilation of antiprotons with baryons is dominated by multi-meson final states 
X. For the parametrisation of the annihilation cross section cr^°'^(s) we take the form given 



in [|T4| for the process 

p + B — B = N,A,... (15) 

Using the same philosophy as for the calculation of the production rate, a simple dif- 
ferential equation for the decrease of the antiproton density in phase space can be written 
down: 

This equation shows that the absorption rate is given by the thermal average of the annihila- 
tion cross section A[^BiT] = {va^^'^) where v is the p velocity. Although a^^'^ls) is strongly 
peaked for low-energy antiprotons, the thermal average {va^^'^) is practically independent 
of the p velocity in the medium. Therefore, in a stationary medium the shape of the p 
spectrum is maintained; the momentum independence of A leads only to a renormalisation 
of the spectrum by annihlation. Any modifications of the spectral shape must be due to the 
dynamical evolution of the hadronic medium. 

The absolute value of A decreases slightly with increasing temperature; it is of the order 
of 0.5 (fm/c)~^ at normal nuclear matter density. This means that 99% of the antiprotons 
are absorbed within 9 fm/c. Of course, due to the rather large value of the absorption rate A 
which is of the order of the size of the antiproton, the classical approach chosen in our model 
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becomes questionable. More reliable results should be based on a quantum field theoretic 
calculation which is beyond the scope of this paper. 



IV. ANTIPROTON SPECTRA FROM AN EXPLODING FIREBALL 
A. A Model for the Heavy Ion Collision 

In order to compare the results of the two previous sections with experimental data 
we connect them through a dynamical model for the heavy ion reaction. In the spirit of 
our thermodynamic approach the so-called hadrochemical model of Montvay and Zimanyi 
is applied for the simulation of the heavy ion collision. In this picture the reaction is 
split into two phases, an ignition and an explosion phase, and particles which have at least 
scattered once are assumed to follow a local Maxwell-Boltzmann distribution. In addition, 
a spherically symmetric geometry is assumed for the explosion phase. The included particle 
species are nucleons, A-resonances, pions and p-mesons. 

As initial condition a Fermi-type density distribution is taken for the nucleons of the 
incoming nuclei, 

P^^^ = 1 + jr\)/. ' (17) 

with the nuclear radius Rq and a surface diffuseness z = 0.55 fm. Because the Fermi motion 
is small compared to the typical bombarding energy of 1 GeV/A, it can be neglected, and 
the initial phase space distribution therefore reads 

fix, p, t) = 6{p,) 6{py) 5{p, - po) p{\x - . (18) 

Here ^(t) describes the path of the center of the nucleus. 

Due to elastic and inelastic scattering hot and dense matter is formed in the central 
region. Its time- dependent chemical composition is governed by the kinetic equations of the 
hadrochemical model given in Ref. |Q. Their particular form needed here can be found 
in fl^. Some of the resulting distributions are visualized in Fig. § for a central ■*°Ca-^°Ca 
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collision with momentum pz = ±700 MeV in the cm. system. One clearly sees that at 
the moment of full overlap of the two nuclei a dense zone with hot nucleons, resonances 
and mesons (not shown) has been formed. Only in the peripheral regions " cold" target and 
projectile nucleons can still be found. On the other hand, chemical equilibrium of the hot 
collision zone is not reached in the short available time before the explosion phase sets in; 
pions and in particular p mesons remain far below their equilibrium abundances Jl^ . 

It is important to note that due to the arguments given in Section || p production is 
strongly suppressed in this initial stage of the reaction. In our simple model we have in fact 
neglected this early p production completely. The ignition phase is only needed to obtain 
the chemical composition of the hot fireball which is expected to be the dominant source for 
the creation of antiprotons. 

For the subsequent expansion of the spherically symmetric fireball into the surrounding 
vacuum analytical solutions can be given if the equation of state of an ideal gas is taken 
as input |^5[. If excited states are included in the model, an exact analytical solution is no 
longer possible. Because a small admixture of resonaces is not expected to fundamentally 
change the dynamics of the system, we can account for their effect in first order by adjusting 
only the thermodynamic parameters of the exploding fireball, but not the expansion velocity 
profiles. 



There is one free parameter in the model a, which controls the density and the 
temperature profiles, respectively. Small values a — represent 5-function like density 
profiles whereas a oo corresponds to a homogeneous density distribution throughout the 
fireball (square well profile). The time-dependent temperature profiles for two representative 
values of a are shown in Fig. ^ for different times t starting at the time tm of full overlap of the 
nuclei. Clearly, a small value of a leads to an unreasonably high temperature (T ~ 200 MeV) 
in the core of the fireball at the beginning of the explosion phase, and should thus be 
considered unphysical. 
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B. Antiproton Spectra from an Exploding Fireball 



Based on the time-dependent chemical composition of this hadrochemical model we can 
calculate the spectrum of the antiprotons created in a heavy ion collision. Let us first con- 
centrate on the influence of the density distribution in the fireball characterized by the shape 
parameter a. Due to different temperature profiles connected with different a values (see 
Fig. ^) the absolute normalization varies substantially when the density profile is changed. 
For more 5-like shapes an extremely hot fireball core is generated, whereas for increasing a 
both density and temperature are more and more diffuse and spread uniformly over a wider 
area. Because of the exponential dependence upon temperature a small, but hot core raises 
the production rate drastically. This fact is illustrated in Fig. ^ for three values of a. Not 
only the total normalization, but also the asymptotic slope of the spectrum is modified due 
to the variation of the core temperature with a as indicated in the Figure. 

Comparing the dotted lines, which give the pure production spectrum, to the solid lines 
representing the asymptotic spectrum at decoupling, the tremendous effect of antiproton 
absorption in heavy ion collisions is obvious. As one intuitively expects absorption is more 
pronounced for low-energetic antiprotons than for the high-energetic ones which have the 
opportunity to escape the high density zone earlier. Therefore, the finally observed spectrum 
is flatter than the original production spectrum. 

Interestingly, while the baryon-baryon and the pp channels are comparable in their con- 
tribution to p production, the pion-baryon channel turned out to be much more effective for 
all reasonable sets of parameters. This fact is indeed remarkable, because here, contrary to 
the discussion in Section |I|, the pions are not in chemical equilibrium; in our hadrochemical 
model the total time of the ignition phase is too short to saturate the pion channel. The 
meson-baryon channel is thus crucial for understanding p spectra. Only by including all 
channels reliable predictions about the antiprotons can be drawn. 

We did not mention so far that in our calculations we followed common practice and 
assumed a finite p formation time of r = 1 fm/c; this means that during this time interval 
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after a p-producing collision the antiproton is assumed to be not yet fully developed and 
thus cannot annihilate. However, there are some (although controversial) experimental 
indications of an extremely long mean free path for antiprotons before final state interactions 



set in |^6[ . We have tested the influence of different values for the formation time r on the p 
spectrum. Fig. |^ shows that this highly phenomenological and poorly established parameter 
has a very strong influence in particular on the absolute normalization of the spectra, i.e. 
the total production yield. In the light of this uncertainty it appears difficult to argue for or 
against the necessity for medium effects on the antiproton production threshold based on a 
comparison between theoretical and experimental total yields only. 

C. Comparison with Experimental Data 

In all the calculations shown above a bombarding energy of 1 GeV/A has been assumed. 
Experimental data are, however, only available at around 2 GeV/A. At these higher energies 
thermalization becomes more questionable ||2^ , and our simple model may be stretching its 



limits. Especially, the temperature in the fireball core becomes extremely high. In order 
to avoid such an unrealistic situation and in recognition of results from kinetic simulations 
iJIJ^ we thus assume that only part of the incoming energy is thermalized - in the following 
a fraction of 70% was taken. 

Fig. H shows calculations for the antiproton spectrum from Na-Na and Ni-Ni collisions 
at a kinetic beam energy of 2 GeV/nucleon. The calculation assumes a p formation time of 
r = 1 fm/c, and takes for the density and temperature profiles the parameter value a = 1 
which corresponds to an upside-down parabola for the density profile. Comparing with the 
GSI data we see our model features too weak a dependence on the size of the collision 
system; the absolute order of magnitude of the antiproton spectrum is, however, correctly 
reproduced by our simple hadrochemical model, without adjusting any other parameters. 
No exotic processes for p production are assumed. As mentioned in the previous subsection, 
the pion-baryon channel is responsible for getting enough antiprotons in our model, without 
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any need for a reduced effective p mass in the hot and dense medium 0,|5[. The existing 
data do not yet allow for a definite conclusion about the shape of the spectrum, and we 
hope that future experiments p8| will provide additional contraints for the model. 



V. CONCLUSIONS 

Heavy ion collisions at typical BEVALAC and SIS energies are far below the pp- 
production threshold. As a consequence, pre-equilibrium antiproton production in such 
collisions is strongly suppressed relative to production from the thermalized medium pro- 
duced in the later stages of the collision. Therefore, p production becomes important only 
when the heavy ion reaction is sufficiently far progressed, in accordance with microscopic 
simulations [Q. By assuming a local Maxwell-Boltzmann distribution for the scattered and 
produced particles forming the medium in the collision zone one maximizes the p production 
rate (see Fig. |1|). If, contrary to the assumptions made in this work, the extreme states in 
phase space described by the tails of the thermal Boltzmann distribution are not populated, 
the antiproton yield could be reduced substantially. 

We also found that the threshold behaviour of the p production cross section is not only 
crucial for the total p yield, but also introduces structures into the initial p spectrum. This 
might give rise to the hope that by measuring the p momentum spectrum one may obtain 
further insight into the p production mechanism. On the other hand we saw here, using a 
Fokker-Plank description for the later evolution of the distribution function fp{p, t) in a hot 
environment, that these structures are largely washed out by subsequent elastic scattering 
of the p with the hadrons in the medium. In addition, the large annihilation rate reduces the 
number of observable antiprotons by roughly two orders of magnitude relative to the initial 
production spectrum; the exact magnitude of the absorption effect was found to depend 
sensitively on the choice of the p formation time r. 

We have also shown that meson (in particular pion) induced production channels con- 
tribute significantly to the final p yield and should thus not be neglected. We were thus able 
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to reproduce the total yield of the measured antiprotons in a simple model for the reaction 
dynamics without including, for example, medium effects on the hadron masses and cross 
sections 

However, we must stress the strong sensitivity of the p yield on various unknown param- 
eters (e.g. the p formation time) and on poorly controlled approximations (e.g. the degree 
of population of extreme corners in phase space by the particles in the collision region), and 
emphasize the rapidly thermalizing effects of elastic final state interactions on the p momen- 
tum spectrum. We conclude that turning subthreshold antiproton production in heavy ion 
collisions into a quantitative probe for medium properties and collective dynamics in hot 
and dense nuclear matter remains a serious challenge. 
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FIG. 1. A(s,t) at different times t calculated from the model of Ref. The starting point 
is a 5-function at s = 5.5 GeV'^. The dashed line is the asymptotic thermal distribution at t=oo), 
corresponding to a temperature T = 133 MeV. 
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FIG. 2. Antiproton production spectrum for different threshold behaviour of the elementary 

production process (x = ^, |, |, | from top to bottom). 
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FIG. 3. Effective temperature Tgs for three Maxwell distributions with initial temperatures 
To = 10 MeV, 50 MeV and 70 MeV, respectively. 
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FIG. 4. Density distributions p(0, 0, z) along the beam axis of target and projectile nucleons 
for a ^^Ca-^^Ca collision, normalized to po = 0.15 fm~^. The solid lines labelled by "incoming 
nuclei" show the two nuclei centered at ±5 fm at time to = 0. The two other solid lines denote the 
cold nuclear remnants at full overlap time t^, centered at about ±3 fm. Also shown are fireball 
nucleons (long-dashed) and A-resonances (short-dashed) at time tm- 
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FIG. 5. Temperature profiles for a = 0.2 and a = 5 at four different times t = (=beginning 
the explosion phase) and t = 3, 6, and 9 fm/c (from top to bottom). 
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FIG. 6. p-spectra for different profile parameters a. The dotted lines mark the initial pro- 
duction spectra. The asymptotic temperatures at an assumed freeze-out density pf = Pq/2 corre- 
sponding to the solid lines are, from top to bottom, 105 MeV, 87 MeV and 64 MeV, respectively. 
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FIG. 7. p-spectrum for different formation times, for a profile parameter a = 1. The dashed 
hne indicates the original production spectrum. 
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